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1. Introduction 

There are reasons to suspect that the vacuum in quantum gravity may determine a preferred 
rest frame at the microscopic level. However, if such a frame exists, it must be very effectively 
concealed from view. Numerous observations severely limit the possibility of Lorentz violating 
physics among the standard model fields [1]. The constraints on Lorentz violation in the gravita- 
tional sector are generally far weaker. 

To allow for gravitational Lorentz violation without abandoning the framework of general 
relativity (GR), the background tensor field(s) breaking the symmetry must be dynamical. Einstein- 
aether theory is of this type. In addition to the spacetime metric tensor field gab it involves a 
dynamical, unit timelike vector field u a . Like the metric, and unlike other classical fields, the unit 
vector cannot vanish anywhere, so it breaks local Lorentz symmetry down to a rotation subgroup. 
It defines a congruence of timelike curves filling all of spacetime, like an omnipresent fluid, and so 
has been dubbed the "aether". This paper aims to provide a review of what has been learned about 
this theory, the status of observational constraints, and some currently open issues. 

A primary motivation for studying Einstien- aether theory — "ae-theory " for short — is the quan- 
tum gravity suspicion mentioned above. A secondary aim is to develop a viable and reasonably 
natural foil against which to compare gravitational observations, in an era when numerous alter- 
nate gravity theories have already been either ruled out or severely constrained. A third source of 
interest is the theoretical laboratory it offers for studying diffeomorphism invariant physics with 
preferred frame effects. 

The action involving metric and aether is highly constrained by diffeomorphism invariance, 
locality, and the unit constraint on u". The only term with no derivatives is the cosmological 
constant, there are no terms with one derivative (other than a total divergence), and there are five 
terms with two derivatives, 

S = ~l^G I ^ (R+K^ n V a u m V b u n ) d 4 x. (1.1) 
Here R is the Ricci scalar, and the tensor K m b is defined by 

K± = c x g ab g mn + c 2 S«5* + c 3 SZ5» + c 4 u a u b g mn , (1.2) 

where the c,- are dimensionless coupling constants. Like general relativity, pure classical ae-theory 

defined by this action is scale-free. (The metric signature is (H ), the speed of light defined 

by the metric g a b is unity, and the aether is taken to be dimensionless.) The term R a bU a u b can be 
expressed as the difference of the C3 and ci terms, up to a total derivative, so is not independent. 
In computations the unit timelike constraint on the aether is usually imposed by adding a Lagrange 
multiplier term X(g a bU a u b — 1) to the action. The covariant derivative operator V a involves deriva- 
tives of the metric through the connection components, and the unit vector is nowhere vanishing, 
hence the terms quadratic in V a u m are quadratic in derivatives of both aether and metric pertur- 
bations, so the metric and aether modes are coupled. The field equations are written out in detail 
in many of the references, beginning with Ref. [2]. It is noteworthy that the aether stress tensor 
includes second derivative terms arising from the variation of the metric in the connections. In all 
phenomenology work to date, it has been assumed that the aether is aligned at large scales with the 
rest frame of the microwave background radiation. 
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The restriction to no more than two derivatives is motivated by the standard precepts of effec- 
tive field theory [3]: higher derivatives would be suppressed by one power of a small length scale 
for each extra derivative. The natural size of the coupling constants c, depends on unknown physics 
at high energies. If the Planck scale is the only relevant scale then the c, are naturally all of order 
unity. If on the other hand there is an additional scale or scales characterizing the Lorentz violating 
physics, then the c, might be naturally smaller and could differ from each other in order of magni- 
tude. Note that although the c\ term in the action is quartic in the aether field, it makes a quadratic 
contribution to the kinetic terms for metric and aether perturbations when expanded around a flat 
background with a constant aether of unit norm. 

Einstein- aether theory is similar to the vector-tensor gravity theories studied long ago by Will 
and Nordvedt [4], but with the crucial difference that the vector field is constrained to have unit 
norm. This constraint eliminates a wrong-sign kinetic term for the length-stretching mode [5], 
hence gives the theory a chance at being viable. An equivalent theory using the tetrad formalism 
was first studied by Gasperini [6], and in the above form it was introduced by Jacobson and Mat- 
tingly [7]. Related vector-tensor gravity theories are that of Kostelecky and Samuel [8] which corre- 
sponds to the Maxwell-like special case of ae-theory (C3 = —c\, C2 = c<\ = 0), both with a fixed norm 
and with a symmetry breaking potential for the vector, of Gripaios [9] which has all two-derivative 
terms and a symmetry breaking potential for the vector, and the generalized Einstein-aether theories 
of Zlosnik, Ferreira and Starkman [10] and of Zhao [11], in which the aether kinetic terms of (1.1) 
are replaced by functions thereof. Kanno and Soda [12] introduced dependence of the coupling 
parameters on a scalar field that has its own dynamics, and studied inflationary cosmology for one 
such model. Other theories involving scalar fields in addition to a vector field with a Maxwell-like 
action for the vector are Bekenstein's TeVeS [13] and Moffat's STVG [14]. 

1.1 Matter couplings 

Lorentz-violating (LV) effects in the matter sector are produced by couplings of the matter to 
the aether. (More complicated Lorentz symmetry breaking patterns require LV extensions of gen- 
eral relativity with other symmetry breaking fields [15].) Such effects have been highly constrained 
by observations [1], so if they exist they are very weak. In this review I will therefore assume that 
matter couples universally to the metric g^. This assumption is motivated by phenomenology, but 
goes against the precepts of effective field theory since the Lorentz violation in the gravitational 
sector would presumably induce Lorentz violating terms in the matter action via loop effects. Su- 
persymmetry could conceivably provide a natural suppression of these terms, as discussed briefly 
in section 9. 

1.2 Metric redefinitions 

When investigating aspects of the theory that do not involve matter, it can be helpful to exploit 
the metric redefinition g' ab = g a y + (£ — \)u a Ub, which "stretches" the metric tensor in the aether 
direction by a positive factor £. (A negative factor would return a Euclidean signature metric.) The 
action (1.1) for (g' ab , u' a = u a / \f^) takes the same form as that for (gab, u a ), with new coefficients 
c\. The relation between the c\ and c,- was worked out in [16], and is conveniently given in terms of 
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certain combinations with simple scaling behavior: 

c'u = c u 

c 'l23 = C c 123 
C13-! = C(ci3-1) 
c[-c' 3 -l = ^-\ Cl -c 3 -\). (1.3) 

Note that in the absence of matter couplings a one -parameter set of Einstein-aether theories is 
actually pure vacuum GR in disguise [17, 16]: the parameters c- all vanish if cu = cm = 2c\ — 
c\ + c\ = 0, and £ > provided Co < 1. 

2. Newtonian and post-Newtonian limits 

In the weak-field, slow-motion limit ae-theory reduces to Newtonian gravity [18], with a value 
of Newton's constant Gn related to the parameter G in the action (1. 1) by 

Gk= ° . , (2.1) 
I-C14/2 

where C14 = c\ + c 4 . (Similar notation is used below for other additive combinations of the c,-.) 
As long as C14 < 2, the Newtonian limit is thus recovered. If cu > 2 gravity is repulsive. This 
suggests the possibility of resolving singularities by anti-gravity if the coupling coefficients c, were 
not constants [6]. Here I discuss only the theory with constant coefficients. 

All parameterized post-Newtonian (PPN) parameters [19] of ae-fheory except the preferred 
frame parameters 0Ji s a agree with those of GR for any choice of the c,. In particular, the Eddington- 
Robertson-Schiff parameters /3 and 7 are both unity [2], the Whitehead parameter ^ vanishes [20], 
and the five energy-momentum conservation parameters 0:3, £1,2,3,4 vanish simply because the the- 
ory is derived from a lagrangian. The parameter (X2 was found for small c; in [21], and the exact 
values of a\ and (X2 were found in [20] : 

_ -8(c 2 3 + Cl c 4 ) 

ai (ci+2c 3 -C4)(2ci+3c2+C3 + c 4 ) 

«2 = ~Z 7^ s U-3) 

2 Ci23(2-Ci 4 ) 

(This particular way of expressing GC2 was given in [22].) In deriving the PPN parameters, the 00- 
component of the metric perturbation and the 0-component of the aether perturbation are expanded 
to 0(2) in the Newtonian potential, while the 0/-components of the metric and i components of the 
aether are truncated at their lowest order, 0(1.5). 1 

Observations currently impose the strong constraints ai < 10~ 4 and ai < 4 x 10~ 7 [19]. Since 
ae-theory has four free parameters c,-, we may set a\ % % exactly to zero by imposing the condi- 
tions [20] 

c 2 = (-2c?-cic 3 + ci)/3 Cl (2.4) 

c 4 = -c\jc\. (2.5) 

'The 0-component of the aether is expanded explicitly only to 0(1) in Ref. [20], but the contribution of the 0(2) 
part is implicitly incorporated by use of the identity (A. 26). 
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With (2.4, 2.5) satisfied, all the PPN parameters of ae-theory are equivalent to those of GR. 

The parameters a\2 can also be set to zero by imposing C13 = cu = 0, but this case is patho- 
logical, as discussed in section 8. Also discussed there is the case C123 = 0, for which (X2 can 
diverge. Another special case occurs when 2c\ — c\ — c\ = 0, for which ai diverges unless also 
c\ + C1C4 = in which case it is indeterminate. Together these imply cu = 2, so Newton's constant 
(2.1) diverges in this case. 



3. Homogeneous isotropic cosmology 

Assuming spatial homogeneity and isotropy, u a necessarily coincides with the 4-velocity of 
the isotropic observers, and automatically satisfies its field equation. The aether stress tensor is 
constructed purely form the spacetime geometry and is identically conserved. It is just a certain 
combination of the Einstein tensor and the stress tensor of a perfect fluid with energy density 
proportional to the inverse square of the scale factor, like the curvature term in the Friedman equa- 
tion [23, 18], 



T. 



a;ther_ C13 + 3C2 



ab 



G a b ~ ^ 3) R(gab+2u a U b ) 



(3.1) 



2 

The latter contribution plays no important cosmological role since the spatial curvature is small, 
while the former renormalizes the gravitational constant appearing in the Friedman equation, yield- 
ing [18] 

(3-2) 



— l + (c I3 + 3c 2 )/2- 

Since G COS mo is not the same as Gn the expansion rate of the universe differs from what would 
have been expected in GR with the same matter content. The ratio is constrained by the observed 
primordial 4 He abundance to satisfy \G cosmo /G^ — 1| < 1/8 [18]. When the PPN parameters 0^2 
are set to zero by (2.4, 2.5), it turns out that G cosmo = Gn, so this nucleosynthesis constraint is 
automatically satisfied [20]. 



4. Linearized wave modes 

When linearized about a flat metric and constant aether, ae-theory possesses five massless 
modes for each wave vector: two spin-2, two spin-1, and one spin-0 mode. The squared speeds of 
these modes relative to the aether rest frame are given by [24] 

spin-2 1/(1 -C13) (4.1) 
spin-1 (2ci-c? + c|)/2ci4(l-ci3) (4.2) 

spin-0 Ci23(2-ci4)/ci4(l-ci 3 )(2+ci3 + 3c2) (4.3) 

The corresponding polarization tensors were found in one gauge in Ref . [24] and in another gauge 
in Ref. [22]. The former gauge choice for the metric and aether perturbations h a b and v" is hot = 
and v\i = 0. In this gauge the nonzero components of the polarization for waves in the 3-direction 
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are given by 

spin-2 h\2i h\\ = —h%2 (4.4) 

spin-1 v 7 , %= [2c 14 c 2 u /(2c 1 -cj + cj)] 1/2 v I (4.5) 

spin-0 v , h m = -2v , h n = h 22 = -cuvq, h 33 = [2c l4 (l +c 2 )/cwi\ v , (4.6) 

where the index / takes only the transverse values 1 , 2. 

The energy of the spin-2 modes is always positive, while for the spin- 1 modes it has the sign of 
(2c i —c\ + c?i)/(\ — C13), and for the spin-0 modes it has the sign of cu(2 — cu) [25, 22]. (These 
reduce to the results of Ref. [26] in the decoupling limit where gravity is turned off.) Under the 
metric redefinition g a }, — > g a }, -+- (£ — l)u a Ub, (1.3) can be used to verify that the squared speeds all 
scale as l/£ as expected, and the signs of the energy are invariant. 

Several constraints are imposed by stability conditions on the linearized modes. The squared 
speeds correspond to (frequency/wavenumber) 2 , so must be non-negative to avoid imaginary fre- 
quency instabilities. They must moreover be greater than or equal to unity, to avoid the existence of 
vacuum Cerenkov radiation by matter [5]. (The existence of ultra high energy cosmic rays, which 
propagate near the metric light cone, requires that the metric and aether modes propagate at no less 
than an extremely small amount below the metric speed of light.) And the mode energy densities 
should be positive, to avoid dynamical instabilities. 

With the 051,2 = conditions (2.4, 2.5) imposed, the squared speeds expressed in terms of 
c± = c\ ±C3 are 



spin-2 1/(1 -c+) (4.7) 
spin-1 (c+ + c- -c + c_)(c + +c_)/4c+c_(l -c+) (4.8) 
spin-0 c+/3c_(l-c+), (4.9) 

the sign of the spin-1 mode energy is that of (c+ +c_ — c+c_)/(l — c+), and c\ 4 = 2c + c_/(c+ + 
c— ). The conditions that the squared spin-2 and spin-0 mode speeds be greater than or equal to 
unity then restrict c± to the region 

< c+ < 1 (4.10) 
< c_ <c+/3(l-c+). (4.11) 



Some algebra shows that the inequalities < c + < 1 and < c_ also ensure that the squared spin- 
1 mode speed (4.8) is also greater than or equal to unity. They also imply < cu < 2, hence 
the condition for attractive gravity mentioned in section 2 need not be separately imposed. Also 
the energies of the spin-0 and spin-1 modes are then positive. That is, when 011,2 = all of the 
requirements for all of the modes are met and gravity is attractive if and only if and c± are restricted 
by the inequalities (4.10, 4.11) [20]. The allowed range is the shaded region in Fig. 1. 

Interestingly, if the mode speeds are instead required to be less than unity (sub-luminal), then 
the spin- 1 and spin-0 energy densities are negative. Hence not only the Cerenkov constraint, but 
also energy positivity (together with 05i 2 = 0) requires mode speeds greater than unity. 
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Figure 1: The shaded region satisfies the stability, Cerenkov, and positive energy constraints (4.10, 4.1 1) on 
c± after having imposed the vanishing of the PPN parameters a,\,2 by the conditions (2.4, 2.5). The dashed 
curve is the locus where the net radiation power from a source with weak self-field matches that of GR. The 
four dotted curves are contours of fixed cu- 

5. Primordial perturbations 

The behavior of primordial cosmological perturbations in ae-theory has been studied in Refs. 
[26, 27]. The coupling to the aether modifies the spin-2 perturbations in two ways: it changes 
their propagation speed, as already discussed above for perturbations of flat space, and it changes 
their amplitude, since that is determined by the "bare" G in the action (1.1) rather than by the 
measured Newton constant Gn- The result is that, given the same Gn, and assuming the PPN 
parameters a,\ 2 vanish, the power in these perturbations differs from that in GR by the factor 
(1 — ci4/2)(l - co) 1 / 2 . When the stability constraints (4.10, 4.11) are satisfied this factor is smaller 
than unity, hence these spin-2 perturbations are even more difficult to detect than in GR. 

Perturbations of the aether can be decomposed into spin-0 and spin-1 modes, both of which 
decay by themselves as the universe expands. In the presence of a scalar inflaton, however, the 
spin-0 aether mode is sourced and no longer decays away. For single scalar field slow-roll inflation, 
it turns out that the power spectrum is generally modified, but it is unchanged when the PPN 
parameters vanish. 

The late time evolution of the primordial perturbations is modified, since anisotropic photon 
and neutrino stresses in the radiation dominated epoch source the spin-1 mode. This leads to 
modified matter and CMB spectra [27]. The effect is rather small however, and is degenerate with 
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matter-galaxy bias and with neutrino masses. 

Primordial perturbations are thus difficult to distinguish in ae-theory and GR. 

6. Radiation damping and strong self-field effects 

Two-body dynamics and radiation in ae-theory have been analyzed for purely weak fields in 
Ref. [22], and including strong self gravity effects in Ref. [28]. In general the radiation includes 
not only spin-2 waves sourced by a varying quadrupole moment, but also spin- 1 waves sourced by 
quadrupole and dipole moments and spin-0 waves sourced by quadrupole, dipole, and monopole 
moments. As seen in section 4, the waves of different spin propagate in general at different speeds 
relative to the aether rest frame. Direct observation of all of these waves is in principle possible 
with existing gravitational wave antennas [19]. With waves of sufficient amplitude from multiple 
sources, and a network of detectors, the aether rest frame and the wave speeds (4.1)-(4.3) and 
polarizations (4.4)-(4.6) could all be measured. 

6.1 Weak fields 

If the fields are weak everywhere (including inside the radiating bodies), and the PPN param- 
eters ai 5 2 vanish, then the dipole source vanishes and the monopole source is uncoupled to the 
fields. In this regime therefore all radiation is sourced by the quadrupole, as in GR. The net power 
radiated in spin-0, spin-1, and spin-2 modes in this case is given by (G^s/ /5)Q i j, where Qjj is 
the quadrupole moment and si = s/[ci\ is a function of the coupling parameters c,- that reduces to 
unity in the case of GR. Agreement with the damping rate of GR (confirmed to ~ 0.1% in binary 
pulsar systems [19]) can be achieved by imposing the condition #/[ci] = 1, which is consistent with 
the constraints (4.10, 4.11), as shown by the dashed line in Fig. 1. 

6.2 Strong field effects 

For compact sources with strong internal fields such as neutron stars or black holes, terms of 
higher order in Gm/d (where m and d are the mass and characteristic size of the source) can play 
an important role in the dynamics and radiation. Since background fields and radiation fields vary 
little over compact source dimensions, these strong field effects can be handled using an "effective 
source" dynamics specified by a worldline action integral [28] 



where v" is the 4-velocity of the body, u a is the local background value of the aether, and o and 
o' are a constants characterizing the body, called a "sensitivity parameters" or just "sensitivities". 
A similar technique was introduced by Eardley [29] for scalar-tensor gravity, in which the first 
sensitivity measures the dependence of the action on the local background value of the scalar field. 
The ellipses in (6.1) stand for all other invariants that can be formed from the particle worldline 
and the background fields. All of these are of higher order in the velocity relative to the aether or 
involve more derivatives and are hence suppressed. Among other things, the sensitivities determine 
the dependence of energy on velocity, E = mo + |(1 + cr)raov 2 + 1(1 + a — o')mov 4 + (This 




(6.1) 
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is obtained by expressing the action (6.1) in the form / dt(T — V) where V = mo, and reading off 
E = T + V.) At lowest order in binding energy Q., a is given by 



where /[c,-] scales as a for a <C l. 2 Note that the leading term vanishes when (X12 vanish. The 
higher order terms in a have not yet been computed, nor have even the lowest order terms in a' 
(which will also scale as a). They presumably differ for black holes and stars, and will depend on 
the equation of state of the star. 

Nonzero sensitivities lead to a number of phenomena that are constrained by observations, 
including dependence of dynamics on the ambient aether vector, modified post-Newtonian two- 
body dynamics, modified quadrupole sourced radiation, and both monopole and dipole sourced 
radiation. When CLx^ = 0, all of these constraints are met provided the sensitivities are less than 
~ 0.001, which will certainly be the case if c, < 0.01 [28]. To be more precise would require 
knowing the actual dependence of the sensitivities on the c\. It is conceivable that the constraints 
will be much weaker. 

To determine the strong field constraints just mentioned, it is necessary to make some assump- 
tion about the speed V of the observed systems with respect to the background aether frame. It 
was argued in Ref. [28] that neglect of terms involving this speed is justified provided V < 10~ 2 , 
which is easily satisfied for any known proper motion relative to the rest frame of the microwave 
background radiation. 

7. Spherically symmetric stars and black holes 

Unlike GR, ae-theory has a spherically symmetric mode, corresponding to radial tilting of the 
aether, and there is a three-parameter family of spherically symmetric static vacuum solutions [31]. 
If asymptotic flatness is imposed and the mass fixed, there remains a one-parameter family [2, 31], 
whereas GR has the unique Schwarzschild solution. As explained below, the evidence suggests 
that only two solutions in this one-parameter family are astrophysically relevant, at least provided 
the couplings are not too large. 

The total energy of an asymptotically flat solution of ae-theory receives a contribution from 
the asymptotic 1/r falloff of the aether [25, 32]. The aether normalization condition locks this 
contribution to that of the metric, so the energy can be expressed in terms of the metric alone. It 
turns out to be just E = ro/G^, where ro is the coefficient of the 1/r term in the g rr inverse metric 
component 3 and Gn is Newton's constant (2.1). That is, it agrees with the Newtonian mass in the 
Newtonian limit, which is inevitable given that the theory has a Newtonian limit. This is how the 
mass of the solutions discussed below is determined. 

2 This corrects an error in version 1 of Ref. [28], where o" is said to scale as cf. (Also the a prefactor C14 in Eqn. 
(70) should be deleted.) As a result of this correction, the likely constraints on c,- are an order of magnitude stronger, as 
stated here [30]. 

3 This is the asymptotic Misner-Sharp energy [33], covariantly defined by ro = lim r ^ 00 (r/2)(l — g ab r a r j,), where r 
is the area radius. 




(6.2) 
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7.1 Static vacuum aether 

The solution outside a static star is the unique vacuum solution for a given mass in which the 
aether is aligned with the Killing vector [31]. This "static aether" vacuum solution depends on the c, 
only through the combination cu, and was found analytically (up to inversion of a transcendental 
equation) [31]. It has a globally timelike Killing vector, and is asymptotically flat. The spatial slices 
have a minimal area 2-sphere, inside of which the 2-spheres flare out to infinite area at a singularity 
on a would-be Killing horizon with vanishing surface gravity. The affine parameter distance to 
the singularity is finite along radial null geodesies. (The minimal two sphere is a property of 
the vacuum solution only; it never occurs in the presence of a matter source.) Although this static 
wormhole indicates the presence of a negative effective energy density in the field equation (i.e. the 
Einstein tensor has this property), all solutions in this family have positive total mass. In Ref. [34] it 
was found that this static aether solution is stable to linear perturbations under the same conditions 
as for flat spacetime, with the exception of the case c\23 = 0. 

7.2 Neutron stars 

The solution inside a fluid star has been found by numerical integration, both for constant 
density [31] and for realistic neutron star equations of state [35]. The maximum masses for neutron 
stars range from about 6 to 15% smaller than in GR when cu = 1, depending on the equation 
of state. The corresponding surface redshifts can be as much as 10% larger than in GR for the 
same mass. Measurements of high gravitational masses or precise surface redshifts thus have the 
potential to yield strong joint constraints on c\n and the equation of state. 

Sufficiently compact stars can have around them an innermost stable circular orbit (ISCO), 
which occurs in the Schwarzscild solution of GR at the radius 6G^M. The ISCO in the static aether 
solution lies at a radius larger by approximately the factor (1 +0.03ci4), with orbital frequency 
lower by the factor (1 — 0.04ci4). These small differences would likely be difficult to measure in 
practice. 

7.3 Black holes 

For black holes the aether cannot be aligned with the Killing vector, since the latter is not 
timelike on and inside the horizon. Instead, the aether is at rest at spatial infinity and flows inward 
at finite radii. The condition of regularity at the spin-0 horizon (where the outgoing spin-0 waves 
propagate at constant radius) selects a unique solution from the one-parameter family of spherical 
stationary solutions for a given mass [36]. This is because the coefficient of the second radial 
derivative term in a field equation vanishes at the spin-0 horizon, which forces the second derivative 
to blow up unless the rest of the terms also vanish. Generically this regularity condition does not 
hold, but the solution can be tuned to satisfy it. When a black hole forms from collapse of matter, 
the spin-0 horizon develops in a nonsingular region of spacetime, where the evolution should be 
regular. This motivated the conjecture that collapse produces a black hole with nonsingular spin- 
horizon, which has been confirmed for some particular examples in numerical simulations of 
collapse of a scalar field [37]. 

The black holes with nonsingular spin-0 horizons are rather close to Schwarzschild outside 
the horizon for a wide range of couplings; for instance, the ISCO radius differs by a factor (1 + 
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0.043ci + 0.061cj), in the case with C3 = C4 = and C2 fixed so that the spin-0 speed is unity [38]. 
(This expansion is accurate at least when c\ < 0.5. No solution with regular spin-0 horizon exists 
in this case when c\ > 0.8.) Inside the horizon the solutions differ more, but like Schwarzschild 
they contain a spacelike singularity. Black hole solutions with singular spin-0 horizons have been 
studied in Ref. [39]. These solutions can differ much more outside the horizon. Quasi-normal 
modes of black holes in se-theory have been investigated in Refs. [40]. 

8. Special values of c. 

It would be interesting if there existed a symmetry of the action restricting the values of the 
coupling parameters c;, but none has been identified. Values for which the PPN parameters ot\ ,2 
vanish evidently lead to more symmetry in the weak field limit, but this does not appear to ex- 
tend to a symmetry of the full theory. A microscopic theory might predict a particular relation or 
hierarchy between the coefficients, but no such theory is in hand. Hence any prior imposition of 
restrictions would be unjustified. Nevertheless, we may ask from a purely phenomenological per- 
spective whether some simplifying restrictions could be compatible with observations. The answer 
is negative. It appears that, unless all of the parameters are much smaller than unity, then in order to 
be compatible with all of the theoretical and observational constraints (leaving aside the radiation 
reaction and strong self-field effects), most likely none of the a, nor C13, cu, or C123 can vanish. 
We now discuss these special cases in some detail. 

• The first case to be examined in detail [8, 7] was C13 = C2 = c\ = 0, i.e. the "Maxwell action" 
(with the unit constraint on the vector). The PPN result for «2 (2-3) is infinite in this case, 
and the spin-0 mode speed is zero. The perturbation series used in the PPN analysis is thus 
evidently not applicable. Independently of that, however, other problems with this case have 
been identified, such as the formation of shock discontinuities [7, 41] and a possibly related 
instability [34]. This case is equivalent to Maxwell theory in a special gauge if it is further 
restricted to the sector in which the Lagrange multiplier field vanishes (see [42, 7, 43, 44] 
and references therein). However, this equivalence holds only if one abandons the notion that 
the vector is itself physical. The shock discontinuities and instabilities of the vector can then 
be regarded as gauge artifacts, but the theory is no longer that of a Lorentz-violating vector 
field. 

• The mode speeds (4.1) are all equal to unity if C13 = C4 = and c% = ci/(l — 2ci). In this 
case the PPN parameters become a\ = — 4c 1 and (X2 = 0, so the observational constraint 
on a\ imposes the stringent condition |ci| < 2.5 x 10 5 , and the other parameters are then 
similarly small. Positivity of the spin-1 and spin-0 mode energies imposes < c\ < 1. With 
small positive c\, the remaining constraints can be satisfied as well. 

• The case c\ = C3 = is evidently not covered by the existing PPN analysis, since a,\ (2.2) 
is indeterminate and the spin-1 mode speed is zero. (If c? = then also the spin-0 speed 
vanishes. If C4 = then the spin-1 speed is indeterminate and the spin-0 speed is infinite, 
but the energy of both modes is zero.) It appears that the post-Newtonian order of some field 
quantities must be modified in this case, which could be interesting to examine ab initio. 
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• Assuming that (i) dip = 0, (ii) the speed constraints (4.10, 4.11) are satisfied, (iii) the spin-1 
mode speed (4.2) is not infinite with finite energy, and (iv) putting aside the case c\ = C3 = 
which is not covered by existing PPN analyses, the only viable case in which any of a, cu, 
C14, or C123 vanishes is the special case C3 = C4 = 2c\ + 3c2 = 0, with 2/3 < c\ < 1. This large 
value of ci is probably inconsistent with the strong field constraints from orbital binaries, but 
as mentioned above those are not yet precisely known because the sensitivity parameters 
have not yet been computed for neutron stars, so this case is not yet ruled out. All the other 
special cases contradict one or more of the conditions (i)-(iv), which can be seen as follows. 
We halt when a case is reduced to a previous case. 

- ci = 0: ai = then implies also C3 = 0. As stated above, the case c\ = C3 = requires 
a new analysis. 

- C13 = 0: The speed constraint (4.1 1) then implies also c_ = 0, so c\ = C3 = 0. 

- C14 = 0: oti = then implies c\ = c\, so either c + = or c_ = 0. In either case, the 
spin-1 speed diverges unless c\ = 0. 

- C123 = 0: The denominator of the second term in 0L2 vanishes (2.3), so either a.2 diverges 
or, if the numerator also vanishes, it is indeterminate. When (X\ = 0, (X2 becomes 

qCi23(2-Ci 4 ) 

so the numerator vanishes in this case only if C13 = 0. (This also implies cu = 0.) 

- C2 = 0: (Z1.2 = implies (2.4), which can be written as 

C2 = -c+(c+ + 3c_)/6ci. (8.2) 

Thus either c + = (covered above) or c + + 3c_ = 0. In the latter case, either c + or c_ 
is negative, violating the speed constraints. 

- C3 =0: OJi = implies c\ = (covered above) or C4 = 0. In the latter case, (2.4) or 
(8.2) then implies that 2c\ + 3c 2 = 0. 

- C4 = 0: Gti = then implies C3 = 0. 
9. Frontiers 

A number of questions remain open regarding phenomenology of a>theory, as well as more 
fundamental issues. In this section some of these are briefly outlined. 

• The most pressing issue is to compute the sensitivity parameters for neutron stars and black 
holes, so that strong field effects on radiation damping rates, waveforms, and post-Newtonian 
dynamics can be evaluated. To determine the sensitivities for a given compact object one 
might compute (presumably numerically) the exact solution for the object moving at some 
velocity with respect to the asymptotic aether, and match e.g. multipole moments of the 
asymptotic fields to those of the corresponding solution derived from the worldline action 
(6.1). It may be adequate to solve for linear perturbations of the exact, static solutions, rather 
than finding the exact moving solutions, but this may not be any more convenient. 
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• Spherical collapse of a star to a black hole would be accompanied by a pulse of spin-0 
radiation, since the aether must adjust from the static configuration (parallel to the Killing 
vector) to the black hole configuration (flowing across the horizon). This radiation may be 
much stronger than that of ordinary gravitational waves arising from asphericity of supernova 
core collapse. Moreover, it could in principle be measured by gravitational wave detectors, 
as seen from the polarization components (4.6). 

• Rapidly rotating black holes, unlike the non-rotating ones, might turn out to be very different 
from the Kerr metrics of GR. If so, this could provide an interesting observational signature 
to detect or constrain. It may be that the easiest way to determine such solutions would be 
by numerical simulation of axisymmetric collapse. 

• Although the linear perturbations all have positive energy for coupling parameters within in 
a certain range, the total, nonlinear energy has not been shown to be always positive in this 
range. (The positive energy theorem of GR does not apply since the aether stress tensor does 
not satisfy the dominant energy condition.) Even a spherically symmetric result in this regard 
would be useful. The question of energy positivity was addressed in Ref. [41], but attention 
was restricted there to the case where the coupling to gravity is neglected, and only some 
special cases were examined. 

• Supersymmetry could potentially allow for a naturally suppressed Lorentz violation in the 
matter sector [45]. To include gravity in such a scheme would require a Lorentz violating 
supergravity. Can ae-theory be supersymmetrized? The argument [46] that supergravity has 
positive energy and therefore so must GR would presumably extend to ae-theory. Therefore, 
supersymmetrizability is likely only possible at most for those parameter ranges for which 
ae-theory has positive energy. 

• Black hole thermodynamics remains a puzzle in ae-theory. It has not been possible to identify 
the entropy and a corresponding "first law" of black hole mechanics, despite the existence of 
an identity relating mass and angular momentum variations to a horizon integral [32]. The 
difficulty is related to the fact that the aether necessarily diverges at the bifurcation surface, 
so the methods usually employed [47] are not applicable. One might expect that the problem 
is related to the fact that, in a theory with multiple limiting speeds, no unique horizon is 
selected. However, even for the cases mentioned in section 8 for which the mode speeds are 
all equal to unity, the entropy has not yet been identified. 

A further important issue is the status of the generalized second law of thermodynamics. It 
has been argued that in a theory with more than one limiting speed, both quantum Hawking 
radiation [48] and classical energy extraction [49] methods can be exploited to violate the 
generalized second law of thermodynamics. This suggests that there may be some deep 
problem with such theories, including ae-theory (with unequal mode speeds). 

• Lorentz- violating matter fields can be consistently coupled to gravity using ae-theory. In this 
way gravitational back-reaction of such fields is well defined. This has been studied in a 
cosmological setting [50], and the issue of renormalization of the corresponding quantum 
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stress tensor has recently been seriously tackled [51]. Counter-terms renormalizing the ae- 
theory action are generally required. Such analyses should enable the quantum back reaction 
to cosmological expansion as well as to Hawking radiation to be studied. 

• Quantization of ae-theory may be instructive to investigate. The dynamical unit timelike 
vector field determines a time flow with respect to which evolution can be referred. This 
"material reference system" could perhaps be employed [52, 53] to address the "problem of 
time" [54, 55] in canonical quantum gravity. 

• Vacuum ae-theory in 1+1 dimensions might provide a useful simplified model for studying 
quantum effects. The classical theory has been investigated [56], where the full solution 
space was mapped out. The solutions include de Sitter and anti-de Sitter spacetimes, with a 
uniformly accelerated aether that is invariant under a two-dimensional subgroup of SO(2,l) 
generated by a boost and a null rotation. (The aether becomes infinitely stretched on certain 
Killing horizons in these solutions.) The only other solution is a peculiar one with singulari- 
ties and no Killing vectors. 

10. Conclusion 

Einstein- aether theory provides a theoretical laboratory in which gravitational effects of pos- 
sible Lorentz violation can be meaningfully studied without abandoning the generally covariant 
framework of general relativity. It is striking that ae-theory can match observations to the degree 
that it can. While of course starting with four free parameters ci 2,3,4 makes this easier, the freedom 
to choose these parameters might have already been used up in setting the PPN parameters to agree 
with those of GR. Instead, only two parameters need be fixed at that stage, leaving c\ and C3 free. 
Then, despite possessing three distinct types of linearized wave modes, all positive energy, stability 
and vacuum Cerenkov requirements on the modes are met within a large common region of the 
(c\,cs) space. Within this same parameter space, the dynamics of the cosmological scale factor 
and perturbations differ little from GR, and non-rotating neutron star and black hole solutions are 
quite close to those of GR (but may be distinguishable with future observations). 

Additional constraints arise from the observed radiation damping rate in binaries. For systems 
with weak self-fields, a constraint of order 10~ 3 would be imposed on one combination of the 
two parameters c\ and C3. Current observations involve pulsars which possess strong self-fields. 
Presuming that the sensitivity parameters for neutron stars turn out to have the expected magnitude, 
such binaries will exhibit effects constraining both the parameters to be less than around 10~ 2 , in 
addition to the stronger radiation damping constraint on one combination. It may turn out that 
observations of multiple systems with different combinations of sensitivities will constrain both 
parameters separately at the 10~ 3 level. 
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